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ABSTRACT 

Both preheating of the intergalactic medium and radiative cooling of low entropy gas 
have been proposed to explain the deviation from self-similarity in the cluster Lx — 7x 
relation and the observed entropy floor in these systems. However, severe overcooling 
of gas in groups is necessary for radiative cooling alone to explain the observations. 
Non-gravitational entropy injection must therefore still be important in these systems. 
We point out that on scales of groups and below, gas heated to the required entropy 
floor cannot cool in a Hubble time, regardless of its subsequent adiabatic compression. 
Preheating therefore shuts off the gas supply to galaxies, and should be an important 
global feedback mechanism for galaxy formation. Constraints on global gas cooling can 
be placed from the joint evolution of the comoving star formation rate and neutral 
gas density. Preheating at high redshift can be ruled out; however the data does not 
rule out passive gas consumption without inflow since z ~ 2. Since for preheated gas 
£cooi > idyn, we speculate that preheating could play a role in determining the Hubble 
sequence: at a given mass scale, high a peaks in the density field collapse early to 
form ellipticals, while low a peaks collapse late and quiescently accrete preheated gas 
to form spirals. The entropy produced by large scale shock-heating of the intergalatic 
medium is significant only at late times, z < 1, and cannot produce these effects. 



1 INTRODUCTION 



Numerical simulations of dark matter halos reveal an al- 
most universal density profile over a wide range of masses 
(Navarro, Frenk & White 1997; Ghigna et al 2000; Klypin 
et al 2001; Power et al 2002). If the gas faithfully traced the 
dark matter, we would expect the bolometric X-ray luminos- 
ity to scale with the cluster virial temperature as Lx oc T 2 
(Kaiser 1986) , while observations indicate significant steep- 
ening of this relation for low-temperature clusters (Helsdon 
& Ponman 2000, and references therein), suggesting that 
some non-gravitational process is affecting the properties of 
the gas. Indeed, observations of the X-ray surface bright- 
ness profiles in poor groups reveal the existence of an en- 



tropy floor, S = T j; 



2/3 



lOOkeVcm (Ponman, Cannon, 



& Navarro 1999), in the gas. Such an entropy floor reduces 
the central gas density and breaks the self-similarity of gas 
density profiles. There are two schools of thought as to how 
this entropy floor could have arisen: (i) as a result of non- 
gravitational heating of order ~ IkeV per particle, arising 
from winds driven either by supernovae or active galactic nu- 
clei (e.g., Kaiser 1991; Ponman, Cannon, & Navarro 1999), 
(ii) as a result of radiative cooling of the low entropy gas, 
leaving behind only the high entropy gas (which cannot cool 
within a Hubble time; e.g., Voit & Bryan 2001). Numerical 
simulations which incorporate either of these schemes can 
reproduce the observed cluster Lx — T relation. 

In this paper, we critique the first scenario, on the ba- 



sis of the disastrous effects the required level of preheating 
would have on subsequent galaxy formation. Groups and 
clusters are the most massive non-linear structures present 
today; to have a significant effect on their gas profile, a large 
amount of energy must be injected (~ IkeV per particle). 
The virial temperature of galaxies is much smaller (for an L* 
galaxy, typically T v i r ~ O.lkeV); galaxies will therefore be 
unable to accrete gas heated to such high entropies. Further- 
more, we shall show that the gas they do manage to accrete 
cannot cool within a Hubble time. The net effect would be 
a sharp downturn in the comoving star formation rate. Mo 
& Mao (2002) have also recently examined the impact of 
preheating on galaxy formation; they find that preheating 
could have a strong impact on the X-ray luminosity and 
morphology of galaxies. However (unlike the present work), 
they do not consider if preheating will suppress subsequent 
star formation. Throughout, we assume a cosmology where 
(n m ,Q.A,tt b h 2 ,h,a sh -i) = (0.3,0.7,0.019,0.7,0.9), consis- 
tent with current observational constraints (e.g. Netterfield 
et al 2002; Freedman et al 2001; Smith et al 2002). 
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THE IMPACT OF PREHEATING ON GAS 
COOLING 



2.1 The overcooling problem 

Recently, it has been pointed out in semi-analytic work 
(Bryan 2000; Voit & Bryan 2001; Voit et al 2002) and 
numerical simulations (Lewis et al 2000; Muanwong et al 
2002; Dave et al 2002) that radiative cooling could be re- 
sponsible for the observed entropy floor in groups and clus- 
ters. Here we point out that: (i) radiative cooling without 
non-gravitational feedback likely results in severe overcool- 
ing of gas at group scales, where the entropy floor is most 
clearly seen, (ii) Depending on its magnitude and redshift of 
injection, an entropy floor due to preheating can efficiently 
regulate overcooling in groups. Thus, some form of preheat- 
ing is still likely to be necessary. This does not preclude an 
important role for radiative cooling: it reduces the otherwise 
severe energetic requirements for preheating, since the low- 
est entropy gas simply cools, rather than being heated up 
to the entropy floor. 

Let us calculate the amount of cooled gas in the ab- 
sence of feedback, as a function of halo mass. Unlike other 
semi-analytic models, we also wish to consider a dynamic 
model in which gas is drawn inward and compressed within 
a cooling flow, rather than a static model in which all gas 
is cooled out to some cooling radius. This is so as to ex- 
plicitly check that the increased gas densities due to inward 
flow do not significantly increase cooling rates and thus the 
cooled mass fraction. We therefore perform ID spherically 
symmetric calculations very similar to that of Tozzi & Nor- 
man (2001), who compute a sequence of adiabats generated 
by the accretion shock and follow their subsequent evolution 
due to radiative cooling; refer to Tozzi & Norman (2001) 
for more details on such models. This is an implicitly La- 
grangian scheme and allows us to directly track the entropy 
evolution of each mass shell. The main difference between 
our treatment and theirs is that we only implement gas cool- 
ing within a static potential, rather than an evolving poten- 
tial; experimentation has found that this makes little differ- 
ence. We also consider the dark matter to impose a fixed 
potential, and ignore the self-gravity of the gas. 

We approach the problem by considering the sequence 
of adiabats through which a gas shell evolves. The initial 
entropy profile of the gas as a function of the mass pro- 
file, K{M) = P(M)/p(M) 7 = k B T/pm p p 2/3 , can be de- 
termined by using the shock jump conditions if the tem- 
perature and density at the accretion shock are known. We 
assume the dark matter distribution follows an NFW pro- 
file (Navarro, Frenk, & White 1997). We assume a perfect 
gas, P — ^ faT, and a polytropic equation of state, 
P = Kp 1 , where 7 = 5/3. Tozzi & Norman (2001) find 
that the average growth of the main progenitor of the 
dark matter halo, found by running Monte Carlo realiza- 
tions based on the extended Press-Schechter formula (e.g., 
Lacey & Cole 1993), can be approximated by the formula: 

m(z) — ( y+§ ) [ s+Alog ( i+^o )] ^ w here A and B depend on 
cosmology, the final mass M Q , and the observed redshift z a 



(we consider z — in this paper). As a function of (Mdm, 
the infall velocity is then found from: 



V i V S A TTT C S 

f = f + AW '^T 



7-1 



(i) 



where ^ = — (R s w i? v i r is the shock radius, and 
i?ta ~ 2.R v i r is the turnaround radius), p c is the external 
preshock density at R = R s , and c s = (t^igmP? -1 ) 1 ^ 2 is 
the sound speed where Kigm is the preshock entropy of the 
IGM gas being accreted. The postshock temperature is then 
given by (Landau & Lifshitz 1957; Cavaliere, Menci & Tozzi 
1998): 



fc B T; = 



pm p v i 



(1 + VT+i) 



7 

+ io e - 



20 (1 + x /T+e)2 



,(2) 



where e = l5k B T e /4:pm p Vi . The postshock density is given 
by Pi = 9Pc, where the shock compression factor g is given 
by (Cavaliere, Menci, & Tozzi 1997): 



9 = 2(1-^ 



+ 



1/2 



(3) 



From this one can compute the postshock adiabat 

ifahock = ksTi/fimppy' 1 . (4) 

We have also computed the entropy profiles assuming 
that the gas traces the NFW profile and is in hydrostatic 
equilibrium, a model used by Voit et al (2002). We have 
found that the entropy profile computed with either model 
is very similar; both are sharply increasing functions of ra- 
dius. We have also found that the entropy profile for different 
cluster masses (scaled to the entropy of the outermost shell) 
is fairly self-similar across different cluster masses. This is 
not surprising, since density and temperature profiles are 
fairly self-similar in these models; they only differ to the 
extent that their central concentration (which depends pri- 
marily on the collapse redshift) varies. We wil l ex amine the 
effects of an entropy floor in these models in §2.2. 

Given an initial entropy profile K(M,t — 0), 
we simultaneously solve the equations for the variables 
K(M,t),p(M,t),R(M,t) (which automatically also yields 
T(M,i) = K(M,t)p(M,t)~<- 1 pm ii /k B ): 



dP 

dM ga! 

dr 
dM Ea 



GMdm 
4nr 4 
1 

47rp gas r 2 



(•>) 



Given the initial entropy profile we have calculated, we can 
compute the initial pressure profile of the gas, if we assume 
the boundary condition p(Rvii) = /b Pnfw {Rvir ) > where 
/b = fib/f2 m is the universal baryon fraction. The pres- 
sure and entropy profile fully specify the initial density and 
temperature profile of the gas. We now allow for radiative 
cooling to operate, and update the entropy profile via: 



ch> (K) - 



Tcoo\(K) ' 



(6) 



where the cooling time is: 
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(7) 



_ 3 ksT p gas 

''"cool — 7T "7 , 

2 An /im p 

and we use a fit given by Tozzi & Norman (2001) to the 
cooling function An of Sutherland & Dopita (1993), which 
includes both free-free and line emission, for a metallicity of 
Z = O.3/?0. We choose our timestep so that it corresponds 
to the cooling time of the third innermost shell; this is suf- 
ficient to obtain convergent results without excessively long 
computation times. As gas cools and drops out at the cen- 
ter, an inward cooling flow develops. After each timestep, 
we recompute hydrostatic equilibrium in equation Q ac- 
cording to the new profile of adiabats K(M), updating the 
density and temperature profile of the gas. Any gas which 
has cooled below T ~ 5 x 10 5 K is removed from the calcula- 
tion and placed in a cold phase at the center of the cluster; 
we track M coo i and Mhot = M gas ,o — M CO ol- We do not con- 
sider distributed mass drop-out in this calculation. When 
recomputing hydrostatic equilibrium, we use the boundary 
conditions: 



0.6 



r(0) = 
r(M ho t) = r ond 



hot ) — 



"vir 



GM DM (r) 



dr 



(8) 

1/(7-1) 



where p v i v — p (r v i r ), K v i r = K(r\\ Y ) are the initial density 
and entropy of the outermost mass shell. The last bound- 
ary condition uses the fact that the outermost shell is com- 
pressed adiabatically during the inward flow, since for this 
shell, t coo \ ^> tn- Note that since we have 3 boundary con- 
ditions for 2 ordinary differential equations, we can solve 
for r cn d, which is an eigenvalue of the problem. Since this 
is a boundary value problem, we use relaxation to solve the 
equations. The results of our calculations agree very well 
with those of Tozzi & Norman (2000) , although they use a 
slightly different set of equations and boundary conditions. 
In our calculation, we assumed the dark matter potential 
was fixed and merely evolved the equations for a Hubble 
time, allowing the gas to cool. We have also experimented 
with a growing dark matter potential; the results do not 
differ significantly. 

In Figure |l|, we show the fraction / coo i of gas which 
cools in a Hubble time in our models without preheating, 
computed according to the convention: 

r _ Afcooi , Q s 

/coo1 ~m^7Tv W 

where M coo i is the mass in the cool phase and M tot is the 
total gravitating mass. We prefer this to the other widely 
used convention: / coo i = M coo i/(M coo i + M hot ), where M hot 
is the gas mass in the hot phase, because the total gas mass 
Mhot + Mcooi in a fixed halo potential decreases with an 
increasing entropy floor in preheating scenarios (a high en- 
tropy floor hinders shallow potential wells from accreting 
gas). We also choose to show halo mass rather than the X- 
ray temperature because for a fixed potential, the emission 
weighted temperature increases with the entropy floor (since 
gas at the cluster center, which dominates the X-ray emis- 
sion, becomes hotter). These conventions allow the clearest 
depictions of the effect of an entropy floor on gas cooling. 



0.4 



0.2 




cool, global 
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Figure 1. The fraction of gas which cools in a Hubble time / coo i 
as a function of halo mass, as computed in our models with- 
out preheating. The points depict the fraction of baryons in the 
cold phase in a sample of clusters from Roussel et al (2000). The 
dashed line shows the best estimate of the global cold baryon 
fraction from the K-band luminosity function of Cole et al (2001). 
Radiative cooling alone without feedback results in severe over- 
cooling, particularly in groups. 



Also shown as points is / coo i inferred from the data of Rous- 
sel et al (2000), who compile a list of total stellar mass and 
total gravitating mass for a sample of X-ray clusters. As in 
Balogh et al (2001) (see their paper for a very lucid dis- 
cussion of the overcooling problem), we multiply the total 
stellar mass by 1.1 to account for gas in atomic and molec- 
ular form; we also show the global cooled gas fraction of 
/cooi ~ 5% inferred from the K-band luminosity function by 
Cole et al (2001). The latter figure is somewhat lower than 
the customarily quoted value of 6% < ,f co id < 17% found by 
Fukugita, Hogan, & Peebles (1998), but much more robust 
as it does not depend on the highly uncertain relative abun- 
dances of galaxies of different morphological types, which 
have very different M/Lb; by contrast M/Lk varies by at 
most a factor of two over different stellar population histo- 
ries (Bell & de Jong 2000). 

We see that while radiative cooling alone may be consis- 
tent with observations in the most massive clusters, severe 
overcooling of the gas would occur on group scales and below 
if radiative cooling alone were to operate without feedback, 
with ~ 50% of the gas entering the cold phase. In our mod- 
els we have computed the entropy profiles and allowed for 
gas inflow and compression in a cooling flow. We also com- 
puted static models which assume that the gas traces the 
dark matter and that all gas cools out to some cooling ra- 
dius r CO oi where t coo i ~ tn (White & Frenk 1991; Voit et al 
2002; Wu & Xue 2002), and obtained very similar results. 
We comment on this robustness in 



2:2 



This overcooling effect is in fact consistent with all nu- 
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merical simulations which claim to explain the observed 
Lx — Tx relation from the effects of radiative cooling 
alone, although this has not been sufficiently emphasized. 
For instance, while Muanwong et al (2002) find a global 
cooled gas fraction of / CO oi ~ 15%, on group scales where 
M v i r ~ few x lO 13 M0fe _1 , the cooled mass fraction is 
/cool ~ 30 - 50% (see their Figure 3). Dave et al (2002) 
find that while / coo i ~ 24% globally in their simulations, for 
halos with <jid ~ 100km s -1 , then / coc j ~ 80%, while for 
halos with am ~ 500km s _ , then / coo i ~ 50%; for these 
extreme levels of gas cooling, consistency with the observed 
Lx — Tx relation can be achieved. The overcooling prob- 
lem is often dealt with in numerical simulations by delib- 
erately restricting the resolution to ~ 10 1 Mq; all simula- 
tions of higher resolution generically suffer from overcooling 
(Balogh et al 2001), as was pointed out by Suginohara & 
Ostriker (1998). Based on the Mulchaey et al (1996) sam- 
ple, Bryan (2000) claims that the cool gas fraction increases 
from cluster to groups scales, and consistency with the radia- 
tive cooling model can be achieved. However, he computed 
/cooi = M coo i/(Af coo i + Afhot), and estimates of the hot gas 
component are likely to be systematically biased downwards 
in groups due to the limited radial extent of the X-ray profile 
(Roussel et al 2000). 

We therefore contend that radiative cooling alone as an 
explanation for the observed scaling relations and entropy 
floor in groups and clusters is untenable. The problem is 
most acute in groups: since most of the gas in groups is at 
comparatively low entropy, a larger fraction of the gas has 
to cool to recover the observed entropy floor. In the next 
section we explore how entropy injection can ameliorate the 
overcooling problem. 

2.2 The Effect of an Entropy Floor on Gas 
Cooling 

Let us now suppose that some form of external preheat- 
ing takes place, in which the 1GM is boosted to some en- 
tropy floor Kigm- What is the effect of this entropy floor 
on gas cooling? We can easily incorporate the effects of a 
pre-existing entropy floor into our models. If the infall ve- 
locity does not exceed the local sound speed, then the gas is 
accreted adiabatically and no shock occurs; the gas entropy 
.Kigm is therefore conserved. If gas infall is supersonic, then 
the gas is shocked to the entropy 7f s hock computed in the 
previous section. Tozzi & Norman (2001) found that the 
transition between the adiabatic accretion and shock heat- 
ing regime is very sharp. Thus to a very good approximation 
we can set K(M) = max(7f s hock(Af), ^"igm); with this new 
entropy profile we can compute density and temperature 
profiles and gas cooling properties as before. Similar to pre- 
vious work (e.g., Tozzi & Norman 2001; Voit et al 2002), we 
find that an entropy floor creates a core in the density pro- 
file, decreasing the central density and increasing the central 
temperature. 

Let us define, as is customary, the entropy parameter 
S = T/n 2 J 3 , and 

Sl0 ° = (lOOkeVcmO ' (10) 




10 14 10 15 0.1 1 



M nM (Mq h-') S |00 (100 keV cm 2 ) 

Figure 2. The effect of an entropy floor Sioo on the cool gas frac- 
tion /cooi as given by equation H. As the entropy floor increases, 
the fraction of gas able to cool in a Hubble time falls drastically. 
We track the inward flow and compression of the gas as it cools 
(dark solid lines); however, ignoring the work done on the gas by 
the gravitational potential and assuming it is static gives very 
similar results (light solid lines). 

For mean molecular weight fi = 0.6, this is related to our 
previously defined entropy parameter in equation (^) via 
Sioo = 6.27^34, where K34 = 10 34 ergcm 2 g" 5/3 . Typ- 
ical values of the observed entropy floor range between 
Sioo ~ 1 (Ponman, Cannon, & Navarro 1999), to 5*100 ~ 4 
(Finoguenov et al 2002). Only values toward the lower end 
of such estimates can be produced by radiative cooing alone 
(Voit & Bryan 2001), whereas values toward the higher end 
are necessary to reproduce the observed Lx — Tx relation 
(Tozzi & Norman 2001). 

In Figure ^| we show how the cool gas fraction / CO ol 
varies with the entropy floor Sioo ■ We see two main features 
in this plot: (i) an increase in the entropy floor drastically de- 
creases the cooled gas fraction. An entropy floor thus has the 
potential to solve the overcooling problem, (ii) The results 
of the calculation where the inward flow and compression of 
the gas as it cools is tracked (dark solid lines) is very similar 
to the case when this is ignored and the gas is assumed to 
be static (light solid lines). The compression of the gas by 
the gravitational potential as it flows inwards only causes a 
mild increase in the final cooled gas fraction / coo i; naively, 
one might have expected that the increased densities would 
resulted in much more rapid cooling. 

To understand these results, let us consider how gas 
cooling times evolve at fixed entropy. If gas with entropy 
£100 evolves adiabatically, the gas temperature depends on 
its density as: 

T = 4.7 x 10 4 5'ioo5 2/3 (l + 2) 2 K, (11) 
where 6 = n c /n c is the overdensity of the gas. The isobaric 
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cooling time is then: 



^cool 



5n e kBT 
2^n c A(T) 



oc S 3/2 (T < 10 7 K) 
oc — (T>110 7 K) 



(12) 



where fi « 0.6 is the mean molecular weight. We have used 
A(T) oc T 1/2 for T > 10 7 K (free-free cooling dominated 
regime) and A(T) oc T~ 1/2 for 10 5 ' 5 T < T < 10 7 K (line 
cooling dominated regime). We can explicitly compute the 
cooling time as a function of density for a fixed adiabat; the 
results are shown in Figure [| We have used a fit (Nath 2002) 
to the Sutherland & Dopita (1993) cooling function, assum- 
ing that Z — O.SZq ; the detailed temperature dependence of 
the cooling function in the line cooling dominated regime is 
not strictly A oc T _1//2 , which is why there is a weak density 
dependence in the cooling time at low density. The most im- 
portant feature to note is that gas above a certain entropy 
floor 5 cr it > 5*100 cannot cool in a Hubble time, regardless 
of its density and/or how it is adiabatically compressed. It 
is thus easy to understand how an entropy floor drastically 
decreases the cool gas fraction / coo i as Sprchcat — > Scrit- We 
can also understand why the compression of gas in an in- 
ward cooling flow has little effect on / CO oi: the cooling rate 
depends on the entropy of the gas, not the density. We show 
this effect in Fig. Q which shows how the gas density, tem- 
perature and cooling time evolve as a gas parcel at fixed 
entropy is brought inwards. Because gas at large radii is 
at high entropy, an inward flow only brings about modest 
compression; the rise in temperature offsets the cooling rate 
increase due to the increased density, and overall the cooling 
time only shows a very modest decrease. 

The fact that the observed 'entropy floor', Spreheat) cor- 
responds to a critical entropy where t coo i > ts is perhaps 
unsurprising: if icooi(S'prchcat) < tn, then we might expect 
radiative cooling to erase all evidence for an entropy floor. 
Indeed, it was the central insight of Voit & Bryan (2001) 
that since t coo \(S CI it) ~ tu, the entropy floor could have 
arisen from a floor in the radiative cooling time, rather than 
an epoch of preheating. However, these statements have very 
interesting consequences for gas cooling in smaller halos as- 
sociated with L < L* galaxies. For the temperature regime 
we are interested in (T < O.lkeV, or 5 < lOOS 3 ^ (1 + z)" 3 ; 
gas hotter than this cannot be accreted by galaxies), the 
cooling time either stays constant or increases with density 
for isentropic gas, rather than decreasing with density. Thus, 
regardless of its final density profile, gas which is accreted 
adiabatically into a galaxy-size halo has a fairly constant 
cooling time which depends only upon its initial entropy. 
Gas which ends up bound to halos at low redshift will typi- 
cally fulfill 5(1 + z) 3 > 10 at all times. Thus none of the gas 
which is heated to the entropy floor will be able to cool in a 
Hubble time in galaxies. Once the entropy floor in the IGM is 
established, the supply of fresh gas for star formation is cut 
off. This is independent of the details of the accretion of gas 
and halo mergers; in general, shocks due to this processes 
can only increase the entropy of the gas, not reduce it. Thus, 
after z pre heat where the IGM is heated up to some adiabat 
Scritj fresh accretion of gas from the IGM is halted. Here- 




-10 



1000 



<5(l + z) 3 



Figure 3. The cooling time as a function of density n = 2 X 
10 _7 <5(1 + z) 3 cm" 3 for gas at fixed entropy Sioo> as defined in 
equation llOl Also shown are the Hubble time at z = and the 
dynamical time t^ yn oc n~ 1 / 2 . The cooling time depends only 
weakly on density; for the range of temperatures relevant for 
galaxy-sized halos, T < O.lkeV or 5(1 + zf < lOOS 3 ^, it typi- 
cal increases with density. In general, gas above the entropy floor 
S > Sioo cannot cool in a Hubble time. 



after, we shall follow through on the logical consequences 
of the observed entropy floor in groups and clusters, if it 
arose from an epoch of preheating, to assess whether such a 
scenario is realistic or not. 



2.3 Bounds on the epoch of preheating 

We can use the value of the entropy floor to obtain a lower 
bound on the epoch of preheating. The gas must be pre- 
heated at a sufficiently high redshift that it is accreted adi- 
abatically: i.e., Sp rc h C at > •S'shock, where S^ock is the en- 
tropy acquired due to the virialization shock in the ab- 
sence of preheating. Otherwise, the gas will be shocked to 
a higher adiabat, erasing the signature of preheating. The 
entropy acquired at the accretion shock S'shock = T/n 2//3 oc 
Tshock(l + z) -2 where T s hock increases with time as the depth 
of the potential well grows. Therefore, without preheating 
the entropy of freshly accreted gas increased monotonically 
with time. There is thus a minimum redshift z prc h C at below 
which preheating leaves no imprint on the entropy profile of 
a cluster, as the entropy of the gas is boosted above the en- 
tropy floor by the strong virialization shock. For a present 
day cluster of mass M , we can compute its entropy pro- 
file if we know its accretion history and mass as a function 
of redshift M(z). We can use this to bound the epoch of 
preheating: if S prc hcat ~ lOOkeVcm 2 , then for the median 
accretion history of groups with T ~ IkeV (where the en- 
tropy floor is most prominent), the epoch of preheating must 
be z > 2, when the group acquired its core. 
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Figure 4. Evolution of the density, temperature and cooling 
time as a parcel of gas is brought inwards adiabatically in a 
AfnM = 10 14 Mq/i — 1 halo. Because gas at large radii is at high 
entropy, an inward flow only brings about modest compression; 
the rise in temperature offsets the cooling rate increase due to 
the increased density, and overall the cooling time only shows a 
modest decrease. In fact, for adiabatic compression the cooling 
time is almost independent of the density, and depends only on 
the entropy. 



A more detailed calculation can be carried out using the 
extended Press-Schechter theory. We compute the merging 
history of several 1O 14 M halos identified at z = 0, using 
the algorithm of Cole et al (2000), and locate at each red- 
shift the most massive progenitor. For each progenitor we 
compute S'shock, and identify those progenitors for which 
Sshock > Spreheat- These progenitors are unaffected by pre- 
heating. Therefore we are able to estimate the fraction of 
clusters whose most massive progenitors are unaffected by 
preheating as a function of redshift. This fraction is shown 
in Figure |fj| At z — 2.2 approximately 50% of cluster pro- 
genitors are affected by preheating, while at z — 3, approx- 
imately 95% of cluster progenitors are affected. Thus, pre- 
heating must occur prior to z = 2-3 in order to have a 
significant effect on the cluster gas. 

Similarly, from their sample of 9 groups of galaxies ob- 
served with ASCA and ROSAT, Finoguenov et al (2002) 
conclude from extended Press-Schechter theory that the 
accreted gas reaches an entropy level of 400 keVcm 2 by 
z ~ 2.0 — 2.5, while such high entropies were not present 
at z > 2.8 — 3.5, implying that preheating happened in a 
fairly instantaneous fashion at z ~ 2 — 3. 

If the entropy of the IGM was indeed boosted to such 
high levels at z ~ 2, this implies that the supply of cold 
gas to galaxies was shut off at that epoch. Remarkably, in 
their models of cosmic chemical evolution, Pei, Fall & Hauser 
(1999) concluded that the data was consistent with no gas 
inflow since z ~ 2. However, their models were computed 



Figure 5. The mass fraction of a 10 14 Mq halo at z = in ha- 
los which are unaffected by a preheating entropy of Sprchcat = 
lOOkeV cm 2 as a function of redshift. In order to affect a signifi- 
cant fraction of the cluster's mass, preheating must occur before 
z PS 1-2. 



for a flat SCDM model, which has since fallen out of favour. 
By performing a stripped down version of their models, we 
can nonetheless ask whether the data are consistent with no 
gas inflow since 2 pre heat- If fresh gas cannot cool, then 



fi g = -fi,(l - R), 



(13) 



where fi g , Q* are the comoving densities of cold interstellar 
gas and stars, all in units of the present critical density, 
and R is the recycled fraction, which is 0.3 — 0.4 for typical 
IMFs. In this case, the comoving gas density at the epoch 
of preheating is: 



Og(z prehea t) = fig(0) + (1-R) 



cheat J j. 

dz — tl*. 

dz 



(14) 



In principle this is only a lower limit since it neglect outflows 
due to supernova explosions, galactic winds, etc. To model 
the evolution of the comoving star formation rate we use a 
fit to the data from Porciani & Madau (2001) for a ACDM 
cosmology: 



-Rsfr(z) = C(z)h 65 



exp(3.4z) 
exp(3.4z) +22" 



-M 



©yr 



_1 Mpc" 



(15) 



where C(z) = Ah 65 (n m {l + z) 3 + « k (l + z) 1 + Qa) v 7(1 + 
z) 3 / 2 is a correction factor for a ACDM cosmology. We use 
this fitting function largely to model the precipitous drop 
in the comoving star formation rate at z ~ 1; we choose 
the normalization factor A so that (1 — R) J f2„ = /*fib = 
2 x 10 -3 , where /* = 5% is the total baryon fraction in 
stars seen today (Cole et al 2001), as inferred from the 
K-band luminosity function. They also present two other 
models which vary the amount of star formation at high 
redshift; however, we shall see that these variations make 



© 0000 RAS, MNRAS 000, 000-000 



Entropy Injection as a Global Feedback Mechanism 7 



little difference to our calculations. To estimate the comov- 
ing density of cold gas, we use estimates of the neutral gas 
in damped Ly-alpha systems (Pefoux et al 2002; Rao & 
Turrishek 2000), as well as estimates of the local HI density 
from 21 cm emission surveys (Zwaan, Briggs & Verheijen 
2001; Rosenberg & Schneider 2002). 

In Figure m, we show the results of computing fig(z) 
from equation (H). This shows the comoving density of in- 
terstellar gas required at redshift z, if all gas inflow/cooling 
ceases at that redshift, in order to be consistent with the 
observed comoving star formation rate as a function of red- 
shift and the local observed gas density. One simple way to 
phrase this constraint is that f2g(z pr eheat) ~ fi g (0) + fi*(0) « 
2 x 10~ 3 (/ C ooi/0.05), where / coo i is the cool global baryon 
fraction today. We see from the figure that there is suffi- 
cient cool gas at high redshift to produce the stars seen 
today; in particular, preheating at high redshift z ~ 2 can 
be plausibly accomodated without subsequent gas cooling 
and inflow. Note that the normalization of Q B (z) is fixed by 
the assumed value of / CO ol(z = 0) (which is uncertain to at 
least a factor of 2); if this value increases all the curves shift 
upward. Also, Odla is obviously merely a lower bound on 
ficool- The surveys are based on optically selected quasars, 
so dusty damped Ly-alpha systems will not be represented 
if they dim the background quasar so that it falls out of a 
magnitude limited survey; also, a signficant fraction of the 
gas could lie in lower column density systems or in molec- 
ular or ionized form. Therefore, despite uncertainties in the 
modelling and data, we can at least conclude that the data 
do not rule out the absence of gas cooling and inflow into 
galaxies since z ~ 2. 



3 DISCUSSION 

3.1 Shock heating of the Warm Hot Intergalactic 
Medium 

Cosmological simulations predict that ~ 30 — 40% of present 
day baryons reside in a warm-hot intergalactic medium 
(WHIM), with temperatures 10 5 < T < 10 7 K, and mean 
overdensities S ~ 10 — 30, rather than in virialized objects 
(Cen & Ostriker 1999; Croft et al 2001; Dave et al 2001). 
Can gas heated in such large scale shocks accrete and cool in 
galaxies? Dave et al (2001) find that the WHIM at z=0 ap- 
proximately follows the equation of state: p/pb = T/W 4 7 K. 
We can use this to estimate the entropy of the WHIM: 



Swhim(z = 0) 



: 340 ( A 
V30 



1/3 



keVcm 



(16) 



Since Swhim > Scot, none of the gas in the WHIM, even 
if accreted onto galaxies, can cool in a Hubble time. Thus, 
shock heating by large scale structure could conceivably play 
a similar role as preheating in raising the entropy of the IGM 
and suppressing accretion and gas cooling. 

However, in hierarchical structure formation scenarios, 
such shock heating cannot be responsible for the observed 
entropy floor in groups and clusters. Shock heating of the 
WHIM corresponds to ~ la fluctuations turning nonlin- 
ear today; the cores of groups and clusters were assembled 
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Figure 6. The comoving density of interstellar gas required at 
redshift z, if all gas inflow/cooling ceases at the redshift, in order 
to be consistent with the observed comoving star formation rate 
as a function of redshift, as computed from equations 
The curves are normalized to produce the total comoving den- 
sity of stars seen at z=0, fi*(0) fa 2 X 10 -3 , as inferred from 
the K-band luminosity function (Cole et al 2001); any increase in 
this estimate would move the curves upward. The points show 
the inferred comoving neutral gas density in damped Lyman- 
alpha systems (Peroux et al 2002; filled circles; open circles denote 
their correction for gas in lower column density sytems), Rao & 
Turnshek 2000 (filled squares) and seen locally in 21 cm emission 
(Zwaan et al 2001 (filled triangle), Rosenberg & Schneider 2002 
(open triangle)). In general these represent lower limits. There is 
sufficient cold gas at high redshift to be consistent with a shutoff 
in gas accretion since z ~ 2. 



at high redshift when they were > 2a fluctuations. Gas 
shocked in the WHIM phase is likely shocked again as it is 
boosted to the higher entropies S ~ 1350(T/lkeV)(n/2 x 
10 _5 cm _3 ) _2 ' /3 keVcm _2 associated with the outskirts of 
present-day groups and clusters. Shock heating of the 
WHIM is likely only to affect gas accretion in voids, where 
the halos have shallow potential wells and are unable to ac- 
crete gas of high entropy. 

We can estimate the redshift evolution of the WHIM 
entropy as follows. The postshock temperature of the WHIM 
can be estimated as: 



Twhim(z) = K(HL n i 



(17) 



where H(z) is the Hubble constant, L n i is the nonlinear 
length scale in proper coordinates; this produces a good fit to 
the results of numerical simulations (Cen & Ostriker 1999; 
Dave et al 2001). We have chosen the normalization con- 
stant K so as to reproduce the median temperature for a 
given overdensity obtained by Dave et al (2001) at z=0. 
At a given overdensity S, the entropy of the WHIM scales 
as Swhim oc H{z)' 2 L 2 ll {l + z)~ 2 \ it thus falls rapidly with 
redshift. We plot Swhim (z) in Figure (R); it only exceeds 



© 0000 RAS, MNRAS 000, 000-000 



8 Oh & Benson 




1 2 



Figure 7. The evolution of the median entropy of the WHIM as 
a function of redshift. Note that SiJwhim °c (5 1 / 3 ; here we assume 
8 ~ 30. Although in practice there will be wide scatter about 
this relation, we nonetheless see that the entropy of the WHIM 
falls very rapidly with redshift. In particular, it only exceeds the 
critical entropy S CI n required to prevent cooling within a Hubble 
time at late times. 



Scrit at z < 1. In practice there will be significant scatter 
about the median entropy of the WHIM (for instance, the 
temperature distribution broadens at high redshift; see Fig. 
4 of Dave et al 2001). It would be interesting to use existing 
numerical simulations to compute the distribution function 
of entropy in the WHIM, as well as the redshift evolution of 
the median IGM entropy; such quantities to date have not 
been calculated. 

It is amusing to speculate on the effect of a cosmological 
constant on the future star formation history of our universe. 
The growth function is given by (Heath 1977): 

^"^f^ (18) 

where a = 1/(1+2) is the scale factor and H(a) = (f2 m a -3 + 
flka~ 2 + Oa) 1,/2 , where fik = 1 — fim — ^a- Note that at 
early times when a is small, Di(a) oc a. If we normalize 
this so that a — 1 and -Di(O) = 1 today, then in the future 
when a — > oo, D\(a) — > 1.39. Thus, structure formation is 
already freezing out; the non-linear mass M» will increase by 
at most a factor of 4 from its present day value, in contrast 
to its rapid growth in the past. Thus, gravitational shock 
heating is slowing down and becoming unimportant, as the 
universe enters a period of exponential expansion. Gas in 
already virialized halos will eventually be able to cool as 
the universe ages, modulo feedback by star formation and 
AGN activity fueled by gas cooling. This is in contrast to the 
fl m = 1 SCDM case, when D\ = a and structure formation 
and gravitational shocks continue apace; in this case most 



of the non-stellar baryons in the universe would end up in a 
hot phase which will never cool. 



3.2 Could Preheating determine Galaxy 
Morphology? 

We speculate that preheating may be responsible for the ori- 
gin of the Hubble sequence. Gas in dark matter halos which 
form before the epoch of preheating are assumed to form 
an elliptical galaxy in a rapid collapse, since halos below 
group scales have gas with cooling times shorter than their 
dynamical time t coo \ < tg (Rees & Ostriker 1977). On the 
other hand, halos which collapse after preheating will ac- 
crete gas with entropies more characteristic of groups and 
clusters; in particular, they can cool their gas only slowly, 
with t coo \ > iff. Cooling gas in such halos may instead form 
a disk galaxy. In fact, Mo & Mao (2002) point out that 
preheating could potentially solve two problems associated 
with disk galaxies: explaining the anomolously low X-ray 
surface brightnesses associated with their halos (since pre- 
heating decreases gas densities and X-ray emissivities), and 
preventing angular momentum transport from the disk to 
the dark matter during disk formation (because the gas dis- 
tribution is more extended than that of the dark matter, it 
retains a higher specific angular momentum during merg- 
ers). We note in passing that the low X-ray luminosities of 
halos associated with disks is perhaps unsurprising: a simple 
extrapolation of the Lx — Tx relation to halo temperatures 
typical for the hosts of disks (~ O.lkeV) gives X-ray lu- 
minosities consistent with observations (few x 10 40 erg s _1 ). 
Thus, whatever physics is suppressing the X-ray luminosities 
of group operates in disk galaxies too. 

In a hierarchical formation scenario, the progenitors of 
massive clusters typically formed earlier than those of lower 
mass systems. Thus, massive clusters contain a larger frac- 
tion of galaxies which formed prior to the epoch of preheat- 
ing; the fraction of elliptical galaxies should increase as the 
mass of the dark matter halo increases. To compute this frac- 
tion of galaxies which formed before preheating (and there- 
fore formed ellipticals) we create merger trees for halos of 
different masses at z = 0. We then identify when 1O 12 M0 
progenitor halos form in these trees. If the progenitor formed 
before preheating (assumed to occur at z — 2 here) the halo 
is flagged as hosting an elliptical galaxy, while halos which 
collapse later are assumed to host a disk galaxy. In Figure ^ 
we show the fraction of galaxies formed in 1O 12 M0 halos 
which are ellipticals as a function of the mass of the halo in 
which they are found at z = 0. This is compared against 
the observational determinations of Balogh et al (2002) 
of the elliptical fraction (shown as circles), with their es- 
timate of the "field" elliptical fraction shown as a dashed 
line. Our simple calculation fits rather well with these data 
(note that the "field" should correspond approximately to 
M* ~ 10 13 h" 1 Mq halos). Since the typical mass of halos 
increases in denser environments we expect this model to 
produce the observed morphology-density relation, at least 
qualitatively. A direct comparison requires a detailed model 
of galaxy formation, since it requires calculating the local 
galaxy density within the model. 
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Figure 8. The fraction of galaxies formed in 1O 12 M0 halos which 
are ellipticals, as a function of the mass of the halo in which they 
are found at z = 0. Points with crrorbars show the determinations 
of Balogh et al (2002). (The X-ray luminosities of clusters in 
Balogh et al (2002) were converted to cluster virial masses using 
the relation of Reiprich & Bohringer (1999).) The dashed line 
shows the "field" elliptical fraction from Balogh et al (2002), 
with dotted lines showing their errorbars. 

This is of course little more than a consistency check, 
rather than proof of principle. The idea that ellipticals col- 
lapsed at much higher redshift than spirals is an old one. It 
was suggested long ago that there may be a direct correspon- 
dence between peak height and galaxy morphology, with el- 
lipticals forming from rare (~ 3<r) density fluctuations, while 
spirals form from more common (~ 2a) ones (Blumenthal et 
al 1984; Evrard 1989); such a scheme can reproduce their 
observed relative abundances fairly well and would also ex- 
plain the increased clustering strength of ellipticals. Empir- 
ically, it is known that ellipticals exhibit remarkably tight 
correlations in various global properties and have very small 
inferred age dispersion which indicate a high formation red- 
shift z ~ 2 — 3, unless their formation was synchronized to 
an implausible degree. The evidence includes the tightness 
of the fundamental plane (van Dokkum et al 1998), and the 
extreme homogeneity of their optical colors (Bower, Lacey 
& Ellis 1992; Ellis et al 1998). However, the physical mo- 
tivation for early formation of ellipticals rather than spirals 
has been rather weak. 

The interesting feature that preheating introduces is a 
physical reason why there should be a division in morphol- 
ogy between early and late collapse of halos of the same 
mass, because of the change in the cooling properties of 
the gas. There may be some critical level of entropy (likely 
^cooi^crit) ~ Wn) where galaxy formation switches from 
one mode to another. To test these speculations, it would 
be necessary to conduct high resolution simulations of a 
single galaxy, to observe how galaxy morphology and disk 



structure changes as the entropy of the accreted gas is in- 
creased. To some extent, we know already from previous sim- 
ulations (based on suggestions that photo-ionization could 
delay cooling until z ~ 1) that feedback processes which 
delay gas cooling result in less angular momentum transfer 
and more disk-like features (Weil, Eke & Efstathiou 1998). 

Preheating may also help alleviate possible problems 
with the dearth of galaxies found in voids (Peebles 2001; 
Mathis & White 2002; Benson et al 2002). Since the dis- 
tribution of halo masses in voids is shifted to lower masses 
relative to the field, preheating will be particularly effec- 
tive in suppressing galaxy formation in these environments. 
However, lower mass halos also tend to collapse earlier, so 
a larger fraction of these halos will form prior to preheat- 
ing. This latter effect seems to dominate. Using dark matter 
halos from the simulations of Benson et al (2002) , we identi- 
fied halos which formed after preheating at z — 2 (assuming 
the distribution of formations redshifts given by Verde et al. 
(2001) and no correlation of formation time with environ- 
ment). Removing these halos from the sample makes only 
minor differences to the void probability function for ex- 
ample. Although these simulations have a rather poor mass 
resolution for the study of void galaxies, we expect the ef- 
fects of preheating to become even less for lower mass halos, 
since an even greater fraction of these should form prior to 
preheating. 

3.3 Conclusions 

In this paper, we stress that radiative cooling alone can- 
not account for the observed entropy floor in clusters and 
the minimum entropy required to produce the observed de- 
viation from self-similarity in the Lx — Tx relation. This 
is only possible with severe overcooling. Thus, some form 
of entropy injection-perhaps through supernovae or AGN 
winds-seems necessary (though radiative cooling still plays 
a very important role in reducing the energetic requirements 
for producing an entropy floor (Voit et al 2002)). The level 
of entropy injection required to produce observable changes 
in the density profiles of deep potential wells of groups and 
clusters must have a drastic effect on smaller halos which 
host galaxies. We show that entropy injection can indeed 
play a very important role in regulating gas accretion and 
cooling in galaxies, since gas heated to the entropy floor can 
never cool in a Hubble time, regardless of the densities it is 
compressed to: for gas temperatures typical of galaxy halos 
the cooling time depends almost exclusively on entropy and 
is relatively independent of density. Thus, after the epoch 
of preheating, the gas supply to galaxies was cut off. This is 
not inconsistent with observations if an epoch of preheating 
at z ~ 2 is assumed. We speculate that at a critical level of 
entropy i C ooi(>Scrit) ~ iff, the mode of gas collapse changes 
from free-fall and fragmenation (which produces ellipticals) 
to quiescent accretion (which produces spirals). 

There are many aspects of the impact of preheating on 
galaxy formation which would be worth pursuing in greater 
detail, particularly with numerical simulations. One would 
be the redshift and spatial dependence of entropy injection- 
we have only considered a uniform level of entropy injec- 
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tion which appears instantaneously at z ~ 2. Naively, we 
would expect entropy injection to be proportional to the in- 
tegrated star formation, which is proportional to the local 
stellar density and local metallicity. Similarly, the spatial 
variation and distribution function of entropy due to shock 
heating of the IGM deserves attention, and the possibility 
that the entropy of the IGM may be responsible for galaxy 
morphology. It may well be that the entropy of the IGM 
regulates the amount of cooling and star-formation at any 
given cosmological epoch, and is the chief mechanism which 
prevents overcooling at high redshift. 
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